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Abstract. We have performed Hartree-Fock calculations of electronic structure of 
< 10 electrons in a quantum dot modeled with a confining Gaussian potential well. 
We discuss the conditions for the stability of N bound electrons in the system. We 
show that the most relevant parameter determining the number of bound electrons 
is VqR^. Such a property arises from widely valid scaling properties of the confining 
potential. Gaussian Quantum dots having A^ = 2, 5 and 8 electrons are particularly 
stable in agreement of Hund rule. The shell structure becomes less and less noticeable 
as the well radius increases. 
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1. Introduction 

Modern semiconductor technology has allowed to fabricate and manipulate electrons 
confined within regions of nanometer size having a plethora of shapes. The possibility 
of tunning the shape, size and number of bound electrons of those nanostructures has 
raised a great deal of interest on the subject of the confined low-dimensional few-electron 
systems known as quantum dots (QDs) [HIS] either for specific applications or exploring 
new fundamental phenomena at a quantum level [31 HI [5]. 

Quantum dots have been also termed "artificial atoms" because their electronic 
structure and properties resemble those of natural atoms [6J. Electrons in quantum 
dots are confined due to potential barriers, in much the same way as electrons in 
atoms are confined due to Coulomb nuclear attraction. Further similarities arise 
because electrons within quantum dots interact with each other through Coulomb 
forces, their energy spectra present both discrete and continuum states giving rise to 
binding and dissociation processes, while transitions between them give rise to emission 
and absorption of radiation. Thus it is justified, to some extent, to apply -mutatis 
mutandis- methods of atomic or molecular physics to the study of these systems of 
electrons acted upon by a given confining potential. A spherically symmetric confining 
potential can provide a model for semiconductor spherical nanocrystals embbeded within 
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an insulator In the literature dealing with calculations of electronic structure and 
properties of QDs, a confining harmonic potential model is the most ubiquitous one. 
This can be justified on the basis of a generalized Kohn's theorem [8] satisfied by such 
a potential and, approximately, verified by experiments. A potential of this type allows 
to explain the electronic shell structure observed, e.g., in small clusters. Nevertheless, 
the infinite height and range of a parabolic potential is clearly unphysical and allows to 
accomodate an infinite number of bound electrons, what precludes the consideration of 
binding and dissociation processes. Other finite-range models have been also advocated; 
e.g., the energy spectra of two- and three-electron systems in a spherical potential 
well of finite depth was obtained variationally [9J, and the unrestricted Hartree-Fock 
(UHF) method was applied to the same model for the calculation of the electronic 
structure of systems having up to 20 electrons [1^. More recently, the UHF calculation 
of ground state, chemical potential and charging energies of electrons in an infinite 



spherical potential well, with and without magnetic field, was reported pj]. However, 
the sharp discontinuity at the QD radius for such a potential well is not completely 
satisfactory from neither a physical nor a computational point of view, and some 
interpolating potentials [12] as well as smoothly varying potentials have been proposed 
[131 [HI [13 [IHl [H] • Among them, the Gaussian potential has received some attention 
and its one- [TH [191 [SQl [21] and two-electron spectra [HI [151 [IS] have been calculated. 
Furthermore, this potential is particularly suitable for an atomic-like treatment such as 
UHF calculations. The Hartree-Fock method is one of the most widely used methods for 
the calculation of electronic structure of atoms and molecules. Its usual implementation 
consists in the expansion of the unknown orbitals as a linear combination in a given 
basis set. The method thus determine self-consistently the energy eigenvalues and the 
orbitals. In atomic problems, the most common basis sets are spherical or Cartesian 
Gaussian functions because they provide closed expressions for the matrix elements 
of the atomic Hamiltonian. Such an advantage holds even more when the confining 
potential is Gaussian itself. Hence, we take in this work an atomic-like approach for the 
calculation of the electronic structure of few-electron Gaussian quantum dots. Such an 
approach shall be particularly advantageous when dealing with, e.g., systems of coupled 
quantum dots, and also for the analisys of mixed systems, like molecules coupled to 
QDs. Work along that line will be published elsewhere. On the other hand, such a 
methodology is the starting point for a systematic improvement in the treatment of 
correlation effects by, e.g. higher-order perturbation theory or multi- reference methods 



2. Theory 

We consider a system of few electrons confined by a potential assumed in the form of a 
spherically symmetric Gaussian potential well of a typical radius R and finite depth — Vq, 
i.e., V{r) = —Vq exp (— r^/2i?^). It approaches a parabolic behavior around its minimum 
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while goes smoothly to zero at infinity. In Ref [H] the parabolic approximation 

V{r) = -K,e-"/2^' ^ -Vo + Vo^^ = -Vo + \rmoh\ (1) 
was considered for comparison, where 

^ = —Bi = ' (2) 

relates the frequency u of the harmonic oscillator corresponding to the Gaussian 
exponent A = 1/2R^. The Schrodinger equation for a system of electrons confined 
within the Gaussian well and interacting through Coulomb potentials is 



i i i<j ■' 



^ = i?^, (3) 



where m* is the effective electron masses, and k is the medium dielectric constant 
K = e/eo- By introducing donor Bohr radius ao = {Km/m*)aB as the unit of length and 
donor Rydberg Rd = {m* /mK^) Ry as the unit of energy, Eq. ([3]) becomes 



ijj = eijj, (4) 



such that the results can be easily transferred between different materials by properly 
changing ao and Rd', for instance, for a GaAs semiconductor QD, one has = Q meV 
and Rd = 10 nm. 

The existence of bound states of Eq. (jlj) depends on the radius and depth of the 
Gaussian potential. If both magnitudes are small enough, even a single electron cannot 
form a stable bound state. Some physical insight on the stability of one- and few-electron 
QDs can be gained from a simple variational estimate. The discussion will also prove to 
provide a useful starting point for the systematic construction of basis sets used in the 
more complicated many-electron calculations reported below. Let Eq^^ be the ground 
state energy of a A^-electron QD. Let us estimate variationally the ground state energy 
of one electron in a Gaussian potential. Taking into account that, around its center, 
the Gaussian potential resembles a parabolic one, we propose a a normalized s-type 
Gaussian trial function of exponent a 

(^,(r) = (2a/7r)3/^exp(-ar2). (5) 

The deeper the well, the more similar its ground state is to that of the harmonic 
oscillator. The expectation value of the one-electron Hamiltonian then becomes 

2a 



E(a)^3a-V„[^—) . (6) 

where the first and second terms are the mean values of the kinetic energy and the 
confining potential, respectively. Minimization with respect to a, i.e. dE{a)/da = 0, 
gives 

(A + 2a)^ - SAVo^a = 0, (7) 
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which has to be solved numerically for given A and Vq. Eq. ([7]) has five roots but, 
typically, the one that minimizes E{a), is a positive number of the order of unit a = Oopt- 
The other roots can be rejected on a physical basis because whether they are complex 
conjugate pairs, or real but negative or very close to zero and, hence, the corresponding 
trial functions do not describe bound states. Thus, within this variational approach, 
the energy of a single-electron QD is 

E^^) =E(«opt). (8) 

On the other hand, the ground state of a two-electron QD corresponds to one electron 
with spin +i and the other with spin — | and interacting with each other through the 
Coulomb interaction Vee = 2/ri2. Therefore, 

^(2) = 2^^') + J, (9) 

where J is the energy shift of the one-electron levels due to the Coulomb interaction 
between two Gaussian charge densities 



|2 1 |,„ /„ m2 r3„ .3„ _ A r^opt 



J = 2 |(p,(ri)|^ — |^,(r2)|Vrirf^r2 = 4. (10) 
J ri2 V 

The condition of stability of the two-electron QD discussed in Refs. [9l [10] and apphed 
in Refs. < E^^\ is equivalent to 

4'^ < -J- (11) 

We shall discuss thoroughly the conditions for stability below when considering the re- 
sults for few-electron QDs. 



In the calculations reported in this work we have considered both one- and 
few-electron systems. The former were used for benchmark purposes and provide 
a systematic way for calibrating the method. A one-electron Gaussian potential of 
Vq = 400i?D and A = 1 has been considered by various authors for the comparison 
of the methods proposed in the literature. We also address such a case as our first 
calculation to assess the accuracy of our methodology. The first issue to be considered 
is the choice of a suitable Cartesian Gaussian basis set {v^^*^}, {i = 1, . . . ,K), where 

= x"?/"zPexp(-air2), (12) 

and i = m + n + p is the angular momentum of the function. Hereafter we use the 
spectroscopic notation s, p, d, . . . for £ = 0,1,2,..., thus 

{Vf^} = (13) 

The one-electron QDs were solved by direct diagonalization of the Hamiltonian, 
while the many-electron systems were treated with the unrestricted Hartree-Fock (UHF) 
method. The calculations for a given potential, i.e., for a given pair (Vo,A), were per- 
formed with the same basis sets, irrespective of the number of electrons. 
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Table 1. Values of the parameters A and Vq for the Gaussian potentials studied. 
The radii of the quantum dots are given by i? = 1/V2A. Oopt are the variationally 
optimized exponents used to construct the basis sets. 



A (a.u.*) 


Vo (Rd) 


R (ad) 


a opt (a 


1.0 


400 


0.707 


9.370 885 


0.5 


50 


1.0 


2.183 148 


0.5 


15 


1.0 


1.047 992 



The prescription for the choice of the basis is based on the previously discussed 
variational energy expression [Eq. IQ], and the optimal exponent Oopt [Eq- ©]• 

The values of Oopt for the three cases considered in this paper are listed in Table 
[U Those exponents were then used to generate the basis functions with higher angular 
momentum, i.e., we choose = = = . . . = ctopt- This procedure can be justified 
as long as the Gaussian potential is similar to the parabolic potential, whose solutions are 
Hermite polynomials multiplied by Gaussian functions with an exponent independent 
of the principal and angular quantum numbers n and i, thus exactly satisfying the 
prescription given above. By including one Gaussian function with the optimized 
exponent in each angular momentum block up to a maximum value L, i.e. < £ < L, 
the low-lying states are reasonably well reproduced as compared to previously reported 
values. In particular, the ground state energy is correct within 10~^-Rd approximately, 
as expected, because of the variational procedure for obtaining the basis exponents. 
Generally, states having i > L are missing or, if obtained, have the larger errors. 
This is also a consequence of the fact that the low-lying states of the parabolic and 
Gaussian potentials are alike, but those of high angular momentum are not. Hence, to 
get a correct description of high-lying excited states, functions with the corresponding 
angular momentum i have to be included. In atomic and molecular calculations, a useful 
procedure for enlarging the basis set has been the so-called even-tempered criterion [23] 
by which a basis set having various functions of the same angular momentum has their 
exponent in the same ratio, i.e., for a given £, the exponents af'\af\af \ . . . are in 
the ratio af~^^^/a^p = const. In the limit of a large number of functions the basis 
set should become complete, independently of the ratio chosen. We take arbitrarily a 
ratio of 2 to enlarge every block of angular momentum in the original basis set, thus 
including exponents af^ = aopt/2*~^ {i = 1, . . . , K) smaller than the optimized one. A 
value of = 4 has been enough to reach converged energies within 10~^i?£) for all the 
calculations reported here. In summary, the basis set used in all calculations consists 
of a set {(p^s \ (pp\ . . . , \ ■ ■ ■}, with i = 1, . . . ,K = 4 and £ < L = 4, having 140 
Cartesian Gaussian functions. 
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Table 2. Bound states energies of the one-electron spherical Gaussian potentials 
(Vb,A) = (400,1), (50,0.5) and (15,0.5) calculated by diagonalization of the 
Hamiltonian matrix and Numerov's integration method. Energies are given in donor 
Rydberg units Rd and Icnths in donor Bohr radius a^i, and dots stand for positive- 
energy states. 



n i 






Vo=AOO Rd 






Vo= 


50 Rd 


Vo= 


15 Rd 


Dia^ 


s;on. 


Numerov 


Ref. 


19] 


Diagon. 


Numerov 


Diagon. 


Numerov 


Is 


-341 


895 


-341.892 


-341.^ 


S952 


-35.958 


-35.958 


-7.762 


-7.762 


Ip 


-304 


463 


-304.463 


-304.4628 


-27.282 


-27.282 


-3.697 


-3.698 


2s 


-269 


644 


-269.640 


-269.6445 


-19.987 


-19.987 


-1.215 


-1.215 


Id 


-268 


110 


-268.111 


-268.1107 


-19.204 


-19.204 


-0.439 


-0.439 


2p 


-235 


446 


-235.450 


-235.4500 


-13.109 


-13.111 






1/ 


-232 


849 


-232.875 


-232.8753 


-11.780 


-11.785 






3s 


-203 


983 


-203.979 


-203.9835 


-7.800 


-7.800 






2d 


-202 


427 


-202.431 


-202.4313 


-7.009 


-7.010 






ig 


-198 


700 


-198.798 


-198.7983 


-5.111 


-5.122 






3p 


-173 


156 


-173.244 


-173.2443 


-3.173 


-3.179 






2/ 


-167 


797 


-170.639 


-170.6393 


-1.864 


-1.876 






4s 


-145 


372 


-145.373 


-145.3779 


-0.598 


-0.600 






3d 


-143 


741 


-143.809 


-143.; 


B091 










'29 


-138 


045 


-140.135 


-140.1351 











3. Results and discussion 

3.1. One-electron quantum dots 

In Table [2] we have listed the energy eigenvalues calculated for the potentials studied. 
The first one, Vq = 400, has been studied previously [HI [19], [20], [21] and provides 
a measure of the precision of our calculations. The first column refers to energies 
calculated by diagonalizing the matrix representation of the Hamiltonian after expansion 
in the basis set described above. The second column corresponds to a numerical solution 
of the Schrodinger equation by using the Numerov's integration algorithm [23] , that we 
have implemented as an independent verification for the basis set expansion method. 
For the first potential, where results from a number of methods are available, the third 
column referes to Ref. [19] , which has probably the most accurate results reported so far. 
This potential has bound states with angular momenta as high as £ = 7, what makes it 
difficult to reproduce with a basis set having functions with i = through i = 4. This is 
particularly apparent for the energies i?2/ and E2g, having the largest errors. This effect 
of the incompleteness of the basis set, can be removed by adding four extra functions 
for each i value and symmetry with 5 < £ < 7, what amounts a total of 360 basis 
functions. The energies calculated with this enlarged basis become -E2/ = — 170.639i?D 
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and = — 140.133i?£), in good agreement with Ref. [12]. Nevertheless, the need of 
such a more demanding calculation should be evaluated taking into account whether 
those states are to be occupied or because there is a perturbation exciting the electrons 
from low-lying states to them. 

The case Vq = 50Rd shows a good agreement between the direct diagonalization 
and numerical integration, implying that both the methodology and the basis set used 
are accurate enough. The scheme of levels obtained are in agreement with Ref. |14j . 
The system Vq = 15Rd has not been treated previously and it provides a system having 
just a few bound states, in which the parabolic approximation could be poor. 

It also can be seen from Table [2] that the energy of the corresponding parabolic 
potential, hu = ^IW^jm = 5.48, 10, and 40 Rd, for Vq = 15, 50 and 400 Rd 
respectively, represent approximately, the energy difference between the ground and 
first-excited state. If the potential were really parabolic, that energy difference would 
be the same for every pair of consecutive states. That is not the case for the Gaussian 
potential, i.e., the higher the pair of states considered lie, the worse the approximation 
becomes, as already noted by other authors [HI [16] . The basis sets optimized for these 
one-electron calculations were hereafter used for the few-electrons systems. 



3.2. Few-electron quantum dots 

The results of Hartree-Fock calculations of the ground state energy of A^-electron systems 
(A^ = 1,...,10), as a function of the QD radius, for Vq = 15Rd (left panel) and 
Vq = bORn (right panel), are shown in Fig. [H This figure has features alike the one 
obtained with a finite depth square well [10]. For comparison purposes, the variational 
one- and two-electron energies E^^^ [Eq. ©] and [Eq. ©] were also plotted (dashed 
lines) . 

The figure shows several critical radii Rc^^ at which a crossover between the i^f ^ 
and Eq^~^^ ground states occurs. Thus, Ri^^ is the minimum radius to have bound 
electrons; for R > Ri^\ the A^-electron QD becomes more stable than the (A^ — 1)- 
electron one. It should be noted that at small enough radii no bound state exists. At 
radii less than = 0.430^ (for Vq = IhRo) and i?F = 0.230^ (for R = 50Rd) only 
a single electron can be bound, but for R > Rf \ E^'^ < E^^K The radii for A^ = 1 and 
N = 2 are quite similar {R^^^ ~ R^?^). So also happens for i^c^"* ~ Ri^^ . . . ~ R^^ and 
Rf^c^R^^'l 

Most of the general features of the stability of the A^-electron system can be 
understood considering the one-electron energies. An electron confined within a well 
of a typical size R will have a momentum of order h/R due to the uncertainty principle; 
its kinetic energy then becomes of order /2mR^, which approaches zero as R goes 
to infinity. On the other hand, its potential energy — Vq exp(— r^/2/2^) approaches the 
bottom of the well — Vq. Then, in the limit of large radius, the energy of the one-electron 
system goes to — Vq and, since also the Coulomb interaction goes to zero, the energy of 
the A^-electron well approaches — A^Vq. Such a trend can be observed in Fig. 1. 
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At a given finite radius, however, the energies Eq^^ are not equally spaced for 
successive A^, because of the shell structure of the energy levels of the potential well. 
This shell structure is even more strikingly revealed in the left panel of Figure 2, where 
the chemical potential fi{N) is depicted as a function of the QD radius. The chemical 
potential represents the affinity of the well for binding an extra electron. It is the 
equivalent of the ionization potential or the electron affinity in atomic physics. It 
is clearly apparent a grouping of the lines corresponding to the number of electrons 
occupying the same one-electron level. All curves decreases as R increases due to the 
fact that when the electrons becomes less confined, the system approaches a classical 
behavior and it is easier to add a new electron. In the right panel of Figure 2, the 
charging energy, defined as the difference between the chemical potential of two systems 
differing in one electron, i.e. -E^/i2r ~ /^(^ + 1) ~ /^(^)) is depicted as a function of 
the number of electrons in the QD for three potential radii, namely, 0.8, 1.5 and 2 an- 
The charging energy gives a measure of the stability of the system, the larger Echar, the 
more stable the system. For N = 2,5 and 8 electrons, the system presents large values 
of Echar, in correspondence with the number of electrons needed for filling or half-filling 
a shell. The peaks heights also diminishes as the QD radius increases. That is, the 
charging energy as a function of A^, tends to be fiat because of the disappearance of the 
shell structure in this classical limit. 

As an illustration of the relation between the QD depth and the critical radius, 
we have depicted, in Figure HI the curves Vq versus l/Rc for a QD of three interacting 
(upper curve) and non interacting electrons (lower curve). They are the locus of the 
points representing the minimum QD radii for a given depth. Both curves are nearly 
straight fines, thus showing that Vo-R^ is approximately constant. We shall show, in the 
following, that the variational energy of the one-electron Gaussian potential also has 
such a property. By defining the dimensionless variable x = 2a/ X, the energy EQ{a) 
[Eq. ([6])] can be written as 

3/2" 



E(a) = X 



\-2C^ " 



2 Vx+ 1 



Xe{C), (14) 



where C = VqR^ and e is the dimensionless energy depending on C. On the other hand, 
the equation determining the optimal variational exponent can be written as 

^ = £ - 3C— - = 0. (15) 
dx 2 (x + 1)5/2 V J 

The equation e{C) = represents the condition to have at least one bound state, which 
clearly only depends on VqR^. This equation and f|T5|) are simultaneously satisfied if 
X = 1/2 and C = 9^3/8, i.e., a = X/A = l/SR^ and VqR^ = 1.95. Hence, for a given 
Vq, the minimal radius for having one bound electron is rP = ^/iMJVq = 1.40/v^. 

This can be compared to the results we would obtain with the truncated parabolic 
approximation of the Gaussian potential, i.e., V{r) = —Vo + muj'^r'^/2 ifr < R, and zero 
otherwise, with uj given by Eq. ([2]). In such a case, the condition for having at least one 
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bound state is that the zero point energy becomes less or equal than zero. That is, at 
the critical radius, 

£;« = -V'„ + |fi.= -Vo + |fiy^=0. (16) 

thus giving VqRI = 9/4 = 2.25, i.e., R^^ = comparable to the relation ob- 

tained above. If there were no interactions between the electrons, the energy for allo- 
cating the second electron would be the same, because both of them would occupy the 
one-electron ground state, although with opposite spins. The Coulomb interaction be- 
tween the electron pair, however, modifies such an energy by an amount J. Nevertheless, 
J is, for small- and medium-size QDs, much smaller than the energy difference between 
the ground and first excited one-electron states. Thus, the next critical radius Rc is 
near to the first one R'^\ On the other hand, the radius r'^^ includes the effects of both 
the electron-electron interaction and the shell structure because the third electron has 
to occupy the first excited one-electron level. We can estimate this radius by using an 
argument similar to the one above. The first excited state of the harmonic oscillator is 
now (5/2)/ia; above the bottom of the well, thus giving VqRI = 25/4 = 6.25. A fit of 
the plot of Figure m of the form Vq = a + bRJ gives Vq = 6.30i?~^ for the non interacting 
case, which compares fairly well to the estimation from the truncated harmonic oscilla- 
tor, thus showing that the non-interacting picture is qualitatively correct. Nevertheless, 
the corresponding fit of the UHF calculations gives a relation Vq = 7.4 + 9.0-R~^'*^ as a 
consequence of the Coulomb and exchange interactions. 

Furthermore, the fact that VqR^ determines the number of bound electrons is a 
property shared by a wealth of types of potentials. Let if be a Hamiltonian of the form 

H = -V^-Vou{p) (17) 

where the potential can be factored as a product of a typical energy —Vq times a function 
u{p) of the dimensionless variable p = r/R. That includes most of the types of confining 
potential used in previous works, e.g., the square well V{r) = —Vo9{l — p), the parabolic 
potential V{r) = tvjj{r/RY with R = ^Jfij^muj and, clearly, the Gaussian potential. 
Expressing the Laplacian in terms of p, we get 

^ = ;^ [-Vp - y^RMp)\ , (18) 

which has energy eigenvalues en(Vo, R) = eniVoR?) / R^ . Neglecting the electron-electron 
interaction, the condition of stability for a A^-electron system is just < 0, where n 
is the quantum number corresponding to the highest occupied state. It is fulfilled for 
certain values of Vo-R^ = C'j, giving rise to stability regions. 

Therefore, when the electron-electron interaction Vee can be neglected, the plot of 
Vo versus 1/-Rc represents a phase diagram constituted by regions delimited by straight 
lines delimiting the various zones where N electrons can form a bound state. The argu- 
ment does not hold when V^e is not negligible, since it scales with R differently to the 
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Table 3. Electronic configuration and z-projection of total spin Ms for the most stable 
UHF configurations of a A^-electron Gaussian potential of Vq — 50Rd, calculated in 
the present work. The spin projections of the orbital are represented by "f and |, for 
+1/2 and -1/2, respectively. Spin indices are omitted for closed shell configurations. 
Hund Rule is fulfilled for all the configurations shown. 



N 

1 V 


IVls 


J_/lcty Li UlilL/ V^Uliiig,U.i cxLlUli 


1 


1/2 


lS| 


2 







3 


1/2 


Is^lpi 


4 


1 




5 


3/2 


Is^lpi 
ls^lp|lp| 


6 


3/2 


7 


1/2 


Is^lp^lp^ 


8 





Is^lpS 


9 


1/2 




10 








kinetic energy term. 

It is interesting to analyze how the electrons occupy the UHF spin-orbitals as they 
are added into the system. Table [3] shows the most stable UHF configurations, for given 
A^, along with the total spin projection M5, for Vq = bORn and R > -Rc^°\ In all the 
calculations performed, with < 10, Hund rule is fulfilled, i.e., electrons in the same 
shell maximize the projection of the total electronic spin, in agreement with previous 
works. As an example, we consider the UHF levels of a Vq = 15Rd well of i? = Ian- 
Figure [3] shows how the scheme of levels changes as the first three electrons are added 
to it. Panel (b) shows the one-electron levels of the potential with a single electron 
occupying the ground state Is. When a second electron is added to the system, the 
pair can form a state with total spin projection M5 = 1 [panel (a)] or Ms = [panel 
(c)]. Panels (a) and (d) shows the states available for spins +| and — | drawn adjacent 
to each other. The electron-electron repulsion shifts and splits the one-electron levels 
through the Coulomb and exchange interactions. In the Ms = configuration, both 
electrons can occupy the same orbital, while the Ms = 1 configuration requires the 
second electron to occupy the higher p-^ orbital, thus giving a higher total energy. The 
effect of adding a third electron to the Ms = configuration, is shown in panel (d). The 
Is orbital shifts upwards and splits into ls| and Isj, while the Ip^ becomes stabilized. 

4. Conclusions 

In summary, we have calculated the electronic structure of electrons confined in a 
Gaussian QD by using the UHF method in its Pople and Nesbet's form. We gave a 
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Figure 1. Ground state energies i?Q of TV-electron Gaussian quantum dots (N < 10) 
as a function of the dot radii, for depths Vq — 15Rd (left) and Vq = 5QRd (right). 
The cross-over of levels occurrs when E^f^^^^ < Eq^^ at a critical radii i?c^\ starting 
from approximately i?c^-' — 0.43aD and i?c^-' = 0.22a£), respectively. 




-45 



0.5 



1 1.5 

R(aD) 



-350 



0.5 



1 

R(aj^) 



systematic procedure for constructing the basis set for the self-consistent calculations 
from the one-electron problem and with a given precision. The orbital energies provide 
an insight of the occupancy of the levels and show the fulfillment of Hund rule in 
agreement with previous results obtained with other models flU{ [TT] . 

Since the Gaussian potential has a finite depth, binding and dissociation processes 
can occur. The criterion for the stability of a iV-electron system was already discussed in 
Ref. m in terms of its energy ^ as compared to E^'^ . The condition < "^^ 
is equivalent, in UHF calculations, to the condition that the highest-occupied molecular 
orbital (HOMO) is bound, i.e., euouo < 0, as a consequence of the Koopman's theorem 
[23] . If the electron-electron interaction can be neglected, the regions of the (Vb,-R) 
plane, having bound electrons, are delimited by the relation Vq-R^ = const. This 
feature has been proven to be a widely valid property, as a consequence of the scaling 
properties of the confining potential. 
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Figure 2. Chemical potential (left) and charging energy (right) of the A/'-electron 
Gaussian potentials considered in this work. The grouping of lines in the chemical 
potential and the peaks in the charging energy show the shell structure. 
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Both the chemical potential and the charging energy show a shell structure with 
peaks of stability at = 2, 5 and 8, corresponding to filled and half-filled shells. Both 

magnitudes decreases as the confinement range increases. The methodology applied 
in this work to single QDs, can straightforwardly be transferred to the study of the 
electronic structure and properties of QDs arrays. It is also amenable for considering 
correlation effects with standard atomic methods. Currently, work is in progress along 
such lines. 
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Figure 4. Critical lines representing the relation between the potential depth Vq 
and l/[i?c^'']^, where is the minimal radius to have three bound electrons in the 
system. The lower curve represents the critical line for the non interacting problem 
and its fit to the straight line Vq = 6.3ii~^. The upper curve represents the critical 
UHF Une. The dashed line corresponds to the fit Vq = 7.4 + Q.OR-'^-^^. 
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